Lattice Dynamics in PbMgi/ 3 Nb 2 /303 
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battice dynamics for five ordered PbMg 1/ / 3 Nb 2 /303 supercells were calculated from first principles 
by the frozen phonon method. Maximal symmetries of all supercells are reduced by structural insta- 
bilities, battice modes corresponding to these instabilities, equilibrium ionic positions, and infrared 
reflectivity spectra were computed for all supercells. Results are compared with our experimental 
data for a chemically disordered PMN single crystal. 



I. INTRODUCTION 

Lead magnesium niobate PbMg 1 / 3 Nb2/303 (PMN) is 
an ABCvperovskite that exhibits relaxor ferroelectric 
(RFE) properties, such as a broad, frequency sensi- 
tive, dielectric permittivity peak at T max near room 
temperature^ It was originally suggested that RFE 
properties are associated with a "diffuse phase tran- 
sition" that reflects a random spatial distribution of 
Curie temperatures^ originating from chemical inhomo- 
geneities; however, the implied chemical segregation of 
Mg 2+ and Nb 5+ ions on perovskite B-sites has not 
been observed. Rather, transmission electron microscopy 
(TEM) indicates chemical short-range order (SRO) such 
that 1:1 ordered domains, 20-50 A across, are embed- 
ded in a disordered matrix^ Thus, the microstructure 
is characterized by fluctuations of a 1:1 chemical order 
parameter on the 20-50 A length scale, and not by com- 
positional fluctuations about Mg:Nb=l:2, as previously 
suggested^ 

The "random site model (RSM)" for the ordered 
domains (SiLSiSiiS has NaCl-type Mg/Nb ordering with 
one B-site occupied by Nb 5+ , and the other by a random 
mixture of 1/3 Nb 5+ +2/3 Mg 2 +. Note that the phrase 
"random-site model" is only used for an idealized model. 
The phrase "1:1 phase" is used for real crystals that pre- 
sumably have some SRO chemical correlations which are 
deviations from randomness on the Mg 2+ -rich sites. 

Burns and Dacolil observed that the refractive in- 
dex of PMN departs significantly from a linear tem- 
perature dependence below a characteristic temperature, 
Tb ~ 630-ftT (Burns temperature). They interpreted 
these data as indicating that polar nanoclusters (PNC) 
condense at T < Tb', and Tb is now generally ac- 



cepted as the crossover between RFE and paraelectric 
(PE) states. Recent TEM studies 12 ! 13 have attempted to 
correlate chemical SRO with PNC. 

Chemical SRO makes PMN inhomogeneous at a 
length-scale that affects vibrational mode activities. 
Mode activities in inhomogeneous media are determined 
by the range of interatomic forces, typically a few lattice 
parameters, and they do not depend on the wavelength 
of the probing radiation. However, the probing radia- 
tion wavelength does determine the volume in which the 
spectroscopic response is averaged. If, as in PMN, the 
inhomogeneities are smaller than the wavelength, an ef- 
fective medium approach is justified. 

Room-temperature IR reflectivity spectra of PMN 
single crystals were first published in the 1970s>iiii^ 
Karamyan reported three polar modes, characteristic of 
cubic perovskites. Burns and Dacol, however, reported 
many more modes, which they explained as "two-mode 
behavior" pi^ Subsequent studies of PMN ceramics^ 
agreed with Karamyan's spectrum, but the resolution 
was poor. 

The interpretation of Raman spectra from PMN is also 
controversial (e.g^ii^). The Raman spectra of PMN and 
PST are similarjiii^ and it was assumed that the main 
features of both reflect Fm3m symmetry (locally for 
PMN) as in the 1:1 phase. In Fm3m, the Raman active 
modes are Ai g , E gi and 2F2 g : A\ g and E g modes are ob- 
servable in parallel (VV) geometry; 2F2 g modes are ob- 
servable in crossed-polarized (VH) geometry. Polarized 
Raman studies and Raman spectral intensities support 
TEM results in suggesting that PMN has 1:1 chemical 
SRO in a disordered matrix. 

The objectives of this study are: (1) to perform a 
comparative supercell stability analyses in PMN, includ- 
ing symmetry-breaking relaxations [i.e. finding ground 
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states (GS) with symmetries that are lower than those 
dictated by chemical ordering]. (2) to compute the lat- 
tice dynamics for the same set of ordered PMN supercells, 
and to compare the results with experimental data, e.g. 
by comparing simulated IR reflectivity spectra with the 
experimental one. The goal is to find which ordered su- 
percell most closely approximates the experimental case 
of local 1:1 order. 



II. INFRARED SPECTRA 

IR and THz measurements were performed on a PMN 
single crystal, down to liquid helium temperatures, where 
all polar modes are well distinguished. A disk of 9 mm 
diameter and 3 mm thickness was cut from a PMN single 
crystal, and polished for specular IR reflectivity measure- 
ments. A time-domain THz spectrometer was used to 
determine the complex dielectric response e*(u) in the 
submillimeter and near-millimeter ranges. This spec- 
trometer uses femtosecond laser pulses to generate THz 
radiation via optical rectification on a ZnTe single crys- 
tal, with an electro-optic sampling detection technique. 
Low-temperature spectra were taken in transmission con- 
figuration from 100-900 GHz, where the thin (100 /im 
thick) plane-parallel plate was semi-transparent. Room 
temperature dielectric response up to 2.5 THz was cal- 
culated from the THz reflectance data^ The unpolar- 
ized near-normal reflectivity spectra were taken with a 
Fourier transform IR (FTIR) spectrometer Brukcr IFS 
113v at 20 < T < 300 K in the spectral range 20-650 
cm ; room-temperature spectra were measured up to 
4000 cm" 1 . 

IR reflectivity spectra were fitted together with com- 
plex dielectric THz spectra using a generalized-oscillator 
model with the factorized form of the complex dielectric 
function: 



the fit. The results are shown in Table [I] In particular, 
the TOl phonon frequency shifts from 88 cm -1 at 20 K 
down to 54 cm -1 at 300 K. Therefore we can call it soft 
mode (SM). 

The low-frequency part [so called central mode (CM)] 
of the 300K spectrum, below polar phonon frequen- 
cies, was fit to an overdamped three parameter oscilla- 
tor model. This fit does not explain the experimental 
low-frequency permittivity which is one order of mag- 
nitude higher, owing to broad dispersion between the 
audio- and microwave-frequency ranges iSSiSiiSS It is well 
known that the distribution of relaxation frequencies in 
this dispersion reaches up to the submillimeter range at 
room temperaturei22i2£ At 20 K, a new heavily damped 
excitation appears near 30 cm -1 (Table QJ. 
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where the dielectric function is related to reflectivity R(f) 
by 



R{y) = 



(2) 



VTOj and VLOj are the transverse optic and longitudinal 
optic (LO) frequency of the j-th mode, respectively; 7tOj 
and jLOj the corresponding damping constants. The 
high-frequency (electronic) permittivity Eqo was obtained 
from the frequency-independent 300K reflectivity above 
the phonon frequencies. The temperature dependence of 
£oo is usually very small and was neglected in our fits. 
The TOl mode was fit to the more accurate FTIR reflec- 
tivity. Below 30 cm -1 , the THz data are more accurate 
than the FTIR data, so it was given greater weight in 



FIG. 1: Temperature dependence of Raman intensity (cor- 
rected for the temperature factor). 



III. RAMAN SPECTRA 

Raman spectra were excited with polarized light from 
a coherent INNOVA 99 Ar+ laser (A = 514.5 nm) and 
analyzed with a Jobin Yvon T64000 spectrometer, that 
was equipped with a charge coupled device (CCD). Po- 
larized Raman spectra were measured in backscattering 
geometry on the (1x2x3 mm 3 ) sample. 

Incident light was focused to a spot size of ~ 3/im 
in diameter, as measured by optical microscopy. Low- 
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TABLE I: Parameters of the polar phonon modes in PMN, from fits of IR and THz spectra at 20 and 300 K. Mode frequencies 
vroi, vlch and dampings 7tck 5 IhOi are in cm" 1 , Ae; is dimensionless, £00=5.83. 





20 K 


300 K 


Nr. 
IN O 


VTOi 




VLOi 






VTOi 


jTOi 


VLOi 


IhOi 






36.2 


44.8 


62.5 


56.2 


129.2 


24.2 


69.2 






539.1 


1 


65.0 


15.2 


67.1 


18.2 


0.6 












(SM)2 


88.2 


22.0 


122.9 


32.7 


15.7 


54.3 


38.3 


102.0 


17.3 


79.6 


3 


245.5 


76.9 


264.1 


49.0 


11.8 


227.7 


85.1 


282.6 


65.0 


17.7 


4 


272.3 


59.4 


382.1 


60.9 


4.1 


287.7 


55.0 


334.2 


45.0 


0.9 


5 


385.1 


41.2 


406.3 


21.4 


0.1 


336.2 


47.5 


397.7 


7.5 


0.13 


6 


434.0 


93.7 


455.3 


22.0 


0.3 


427.9 


55.7 


442.3 


35.4 


0.2 


7 


547.9 


80.6 


703.0 


22.9 


1.8 


545.4 


105.1 


703.0 


62.5 


2.0 



and high-tcmperature micro- Raman measurements were 
made with a Linkam FDCS 196 cryostat and Linkam 
TS1500 hot stage, respectively. Measured Raman spectra 
were corrected for the Bose-Einstein temperature factor. 

Fig. shows the temperature dependence of paral- 
lel (VV) and crossed (VH) polarized Raman spectra in 
the temperature interval 77-700 K. The main feature is 
that the frequencies of all observed peaks are only weakly 
temperature-dependent. Below 270 K, a partial depolar- 
ization occurs, due to the leakage of intense peaks from 
the VV spectrum to VH. 

All peaks are best resolved at low temperature, though 
most of them can be traced up to 1000 K. As shown in 
Fig. ^the VV spectrum exhibits: a band at 54 cm -1 ; an 
envelope of strongly overlapping bands between 100 and 
350 cm -1 ; a very weak broad band at about 430 cm" 1 , 
that is clearly seen only at low temperature; strongly 
overlapping bands in the 500-600 cm -1 range; and a 
very intense band at 780 cm" 1 . In the VH spectrum the 
low- frequency band is a doublet (45 and 62 cm -1 ) and 
intermediate bands in the 100-350 range cm -1 are also 
strongly overlapping. The intensity of the 780 cm" 1 band 
increases abruptly at about 270 K, just below T max . The 
780 cm" 1 band was interpreted as a fully symmetrical 
(Ai g ) stretching vibration of oxygen octahedra that origi- 
nates from 1:1 chemically ordered regions^ The 500-600 
cm" 1 wide band, usually associated with E g Raman 
active vibration, increases splitting when T decreases^ 
below 350 K, and splits^i in I4mm symmetry but not 
in R3m. Actually, any monoclinic distortions split this 
mode either. Perfect 1:1 chemical ordering implies four 
Raman-active modes: Ai g +E g +2F2 g . The measured Ra- 
man spectra are much more complicated, however, which 
suggests that both chemical SRO, and local symmetry re- 
ductions from lattice instabilities, activate Raman modes 
that are inactive under Fm3m selection rules. 



IV. COMPUTATIONAL METHODS 

Density functional theory (DFT) calculations for PMN 
were done with the Vienna ab initio simulation pack- 



age (VASP)i2£i2& A plane wave basis set was used for 
electronic wavefunctions and ultrasoft pscudopotcntials 
were employed^ Exchange and correlation energies were 
calculated with the local density approximation (LDA). 
Lattice-dynamical force constants were calculated via the 
frozen phonon method. Berry's phase analyses were used 
to compute dynamical charges, as implemented in VASP 
by M. Marsman. 



A. Ordered Supercells 

Ordered supercells with PMN stoichiometry require 
15n atoms (n = 1, 2...). The calculations presented here 
are for three 15 atom cells, and two 30 atom cells (Fig.[5J: 
a) [001} NNM ; b) [IW]nnm; c) [U1]nnm; d) [001]atcc; 
c) [111]jvt- Here: subscript N = a Nb-layer parallel to 
(h, k, I); M = a Mg-layer parallel to (h, k,l); C = a chess- 
board ordered layer, i.e. a Mg^Nb^-layer or (001)i : i- 
layer; C = a chessboard ordered (001)i : i -layer that is 
displaced by [1/2,0,0] relative to an adjacent C-layer; T 
= a (lll)2:i-layer of composition Mg 2 /3Nbi/3, which is 
triangularly-ordered (Fig. The [001] jvee structure 
(d) has (OOl)jvi, layers that alternate with CC double 
layers (NaCl-type blocks). This structure was proposed 28 
as a possible cation ordering GS for PMN, however 
the stability analyses reported below indicate that the 
[111] a/t structure is lower in energy. The [111]jvt struc- 
ture has (lll)jvt layers alternating with (lll)A/g 2/3 Ar6 1/3 
layers, and, as such, is an ordered approximate of the 
1:1 "random-site model". Although labeled "NT" for 
Nb layer/triangular layer alternation, the cation arrange- 
ment in the 2:1 layers depends on the layer's orientation. 
In Fig. [2k, there are triangularly-ordered layers perpen- 
dicular to [111] and [TTl] and striped layers (Fig. 3b) 
perpendicular to [111] and [Til]. 

For each ordered supercell, two space group symme- 
tries, and two corresponding energies, are given in Ta- 
ble [H] Gchem indicates the space group symmetry dic- 
tated by chemical ordering without additional displacive 
instabilities; Gsss is space group of the lower symmetry 
"structurally stable state." The SSS was determined by 
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TABLE II: Supercell space group symmetries G, formation energies AE (kj/mole, where mole is the Avogadro's number of a 
primitive ABO3 unit cell, and relaxed cell parameters.) 



System 


Gchem 


AEchem 


a 


b 


c 


a 


P 


7 


Gsss 


AEsss 


a 


b 


c 


a 


P 


7 


[001]nnm 


P4/mmm 


16.16 


3.977 


3.977 


12.277 


90 


90 


90 


Pm 


12.27 


4.050 


3.965 


12.250 


90 


90 


89.84 


[110] JVJVM 


Amm2 


14.33 


5.664 


17.087 


4.008 


90 


90 


90 


Cm 


5.50 


5.660 


16.998 


4.092 


90 


90.68 


90 


[111] AW A/ 


P3ml 


14.11 


5.650 


5.650 


6.941 


90 


90 


120 


PI 


5.33 


5.682 


5.759 


6.964 


89.97 


90.48 


120.43 


[001] NCC ' 


P4:/nmm 


11.51 


5.646 


5.646 


12.044 


90 


90 


90 


PI 


2.99 


5.645 


5.665 


12.197 


90.15 


89.93 


89.85 


[HI] ATT 


Immm 


7.66 


5.615 


17.050 


8.022 


90 


90 


90 


PI 


0.00 


(5.613 


16.995 


8.085 


89.99 


89.85 


90.00)* 



*Not the primitive PI cell; 60-atom cell allows comparison with Immm lattice parameters. 



TABLE III: Ionic coordinates (in A) in 15 ion ordered super- 
cells of PMN 





[OOIJawm 


[110] JVJVM 


[111] JVJVM 


type 


X 


V 


z 


X 


y 


z 


X 


y 




z 


Mg 


0.00 


-0.03 


0.00 


3.97 


3.97 


0.01 


-0.12 


-0.04 





03 


Nbi 


0.06 


-0.03 


4.03 


-0.01 


4.04 


-0 


10 


-0.05 


-0.10 


3 


99 


Nb 2 


0.06 


-0.03 


8.22 


4.04 


-0.01 


-0 


10 


-0.16 


4.09 


4 


07 


Oi 


-0.14 


-0.10 


2.16 


1.92 


4.16 


-3 


98 


-0.26 


1.86 


4 


04 


2 


-0.23 


-0.09 


6.12 


0.07 


2.07 


-3 


94 


-0.30 


-0.20 


6 


02 


3 


-0.14 


-0.10 


10.09 


6.00 


-0.16 


-3 


92 


1.79 


3.88 


3 


99 


4 


1.99 


-0.12 


0.00 


-0.16 


6.00 


-3 


92 


-0.29 


-0.18 


2 


06 


5 


-0.14 


1.93 


0.00 


4.16 


1.92 


-3 


98 


-0.37 


-2.08 


4 


01 


o 6 


1.94 


0.01 


4.18 


2.07 


0.07 


-3 


94 


1.83 


-0.17 


3 


98 


7 


-0.17 


1.91 


4.23 


0.11 


0.11 


-2 


02 


-0.23 


3.84 


5 


99 


8 


1.94 


0.01 


8.07 


0.01 


4.01 


-1 


99 


-0.23 


5.96 


4 


07 


9 


-0.17 


1.91 


8.02 


4.01 


0.01 


-1 


99 


-2.14 


3.76 


4 


04 


Pbi 


2.17 


2.04 


1.60 


1.77 


1.77 


-2 


27 


2.17 


2.13 


2 


15 


Pb 2 


2.26 


1.98 


6.12 


6.24 


-2.33 


-2 


23 


-1.97 


2.14 


2 


04 


Pb 3 


2.17 


2.04 


10.65 


-2.33 


6.24 


-2 


23 


-2.01 


2.22 


5 


92 



TABLE IV: Basic vectors for the computed structures 



structure 


multiplier 


x y z 


[001]jvjvm 


4 


1.012460 0.001455 0.000000 
0.001417 0.991364 0.000000 
0.000000 0.000000 3.062478 


[011] JVJVM 


4 


-1.002116 2.002737 0.004304 
2.002737 -1.002116 0.004304 
-0.004184 -0.004184 -1.022970 


[111] JVJVM 


5 


-0.004329 -0.809859 0.797156 
0.814385 0.814646 0.000288 
-0.806681 0.807082 0.798775 


[001]jvcc< 


4 


0.999359 -0.996947 0.000397 
1.002589 1.000282 -0.002638 
0.000558 -0.003683 3.049203 


[111] AT 


4 


1.997368 -1.006931 -1.009787 
-1.006824 1.997432 -1.009407 
0.992322 0.992308 0.001922 




c d 




e 



FIG. 2: PMN supercells: [001]jvjvm (a), [110]jvjvm (b), 
[111]jvjvm (c), [001]jvcc (d), [111]jvt (e) 



applying random perturbations to the atomic positions, 
and then relaxing the system until the absolute values of 
all force constants are less than 0.001 eV/A 2 . The differ- 
ence AEsss — lAEchem is referred to as the relaxation 
energy. 

In general, G S ss C G C hem C Pra3m, and AE S ss < 
AEchem- Randomness of the initial perturbation does 
not guarantee that the SSS is in fact the ground-state of 
the system(G sss = Ggs)> but it does guarantee that it 
is, at least, metastable. 

We show here more tables (see Tables ITTT1 to IXIl|) than 
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TABLE V: The dynamical charges of ions in the 15-ion supercells of PMN 





[001]jViVAf 


[110] NNM 


[lll]jVJVM 


ion i 


7* 


7* 

^ixx 


7* 

^lyy 


7* 

^izz 


7* 

^ixx 


^ixy 


7* 

^iyx 


7* 

^lyy 


7* 

^izz 


7* 

^izx 


7* 

^izy 


7* 

^ixz 


7* 

^ixx 


^ixy 




7* 

^iyx 


7* 

^xyy 


Mg 


2.56 


1.90 


1.90 


1.94 


2.60 


-0.09 


-0.09 


2.60 


2.72 


-0.03 


0.02 


-0.02 


2.72 


-0.03 


0.03 


-0.02 


2.72 


Nbi 


7.64 


8.76 


8.76 


8.76 


8.27 


0.19 


0.19 


8.27 


6.93 


0.17 


-0.13 


0.13 


6.93 


0.17 


-0.17 


0.13 


6.93 


Nb 2 


7.64 


8.75 


8.75 


8.76 


8.27 


0.19 


0.19 


8.27 


6.91 


0.19 


-0.14 


0.14 


6.91 


0.19 


-0.19 


0.14 


6.91 


Oi 


-4.62 


-2.15 


-2.15 


-3.04 


-4.15 


-0.08 


-0.29 


-2.60 


-2.33 


-0.06 


-0.12 


-0.08 


-2.34 


0.17 


0.23 


-0.16 


-5.79 


2 


-5.18 


-2.75 


-2.75 


-3.04 


-2.58 


-0.29 


-0.11 


-4.15 


-5.79 


-0.23 


0.16 


0.12 


-2.33 


-0.06 


-0.17 


-0.08 


-2.34 


O3 


-4.62 


-2.15 


-2.15 


-1.97 


-7.05 


-0.50 


0.16 


-2.01 


-2.34 


0.17 


0.08 


0.16 


-5.79 


-0.23 


0.06 


0.12 


-2.33 


4 


-2.91 


-1.76 


-3.99 


-1.97 


-2.01 


0.16 


-0.50 


-7.05 


-3.88 


0.07 


-0.13 


-0.17 


-2.54 


0.04 


0.12 


0.02 


-2.56 


O5 


-2.91 


-3.99 


-1.76 


-3.04 


-2.60 


-0.29 


-0.08 


-4.15 


-2.54 


0.04 


0.17 


0.02 


-2.56 


-0.12 


-0.07 


0.13 


-3.88 


o 6 


-1.94 


-6.84 


-2.18 


-3.04 


-4.15 


-0.11 


-0.29 


-2.58 


-2.56 


-0.12 


-0.02 


0.13 


-3.88 


0.07 


-0.04 


-0.17 


-2.54 


o 7 


-1.94 


-2.18 


-6.84 


-1.78 


-3.59 


0.55 


0.55 


-3.59 


-3.86 


0.06 


-0.12 


-0.15 


-2.54 


0.05 


0.11 


0.03 


-2.56 


o 8 


-1.94 


-6.84 


-2.18 


-6.74 


-2.05 


-0.05 


-0.05 


-2.05 


-2.54 


0.05 


0.15 


0.03 


-2.56 


-0.11 


-0.06 


0.12 


-3.86 


o 9 


-1.94 


-2.18 


-6.84 


-6.74 


-2.05 


-0.05 


-0.05 


-2.05 


-2.56 


-0.11 


-0.03 


0.12 


-3.86 


0.06 


-0.05 


-0.15 


-2.54 


Pbi 


3.43 


3.80 


3.80 


3.96 


3.97 


-0.27 


-0.27 


3.97 


3.85 


-0.18 


0.12 


-0.12 


3.85 


-0.18 


0.18 


-0.12 


3.85 


Pb 2 


3.33 


3.90 


3.90 


3.97 


3.56 


0.32 


0.32 


3.56 


3.79 


0.41 


-0.51 


0.51 


3.79 


0.41 


-0.41 


0.51 


3.79 


Pb 3 


3.40 


3.80 


3.80 


3.97 


3.56 


0.32 


0.32 


3.56 


3.88 


-0.16 


0.10 


-0.10 


3.88 


-0.16 


0.16 


-0.10 


3.88 



TABLE VI: Diagonal elements of the dynamical matrix (in cm 1 ) for structurally stable [001] awm, [110] awm and 
[lll]iVjvAf structures 



structure 


axis 


Mg 


Nbi 


Nb 2 


Oi 


2 


3 


4 


5 


o 6 


o 7 


o 8 


9 


Pbi 


Pb 2 


Pb 3 


[001] NNM 


z 


268 


332 


339 


725 


585 


724 


335 


343 


296 


316 


255 


322 


120 


106 


120 




X 


379 


224 


224 


254 


285 


251 


492 


300 


597 


264 


601 


263 


62 


80 


62 




y 


410 


145 


145 


234 


288 


233 


274 


533 


230 


525 


230 


526 


50 


54 


49 


[110] jVJVM 


z 


347 


233 


233 


291 


293 


253 


253 


291 


293 


465 


604 


604 


67 


72 


72 




X 


329 


283 


286 


613 


310 


654 


263 


296 


602 


314 


263 


262 


97 


84 


87 




y 


329 


286 


283 


296 


602 


263 


654 


613 


310 


314 


262 


263 


97 


87 


84 


[111] AT JVM 


z 


358 


280 


280 


249 


602 


249 


646 


277 


240 


648 


239 


279 


73 


67 


72 




X 


331 


314 


283 


270 


277 


605 


264 


274 


722 


275 


247 


584 


84 


78 


77 




y 


330 


282 


314 


605 


275 


268 


277 


583 


251 


265 


723 


277 


78 


78 


84 



in our submitted paper, in order to save place in the 
journal. 

B. Short-Range Order and Lattice Averaged Spin 
Products 

Chemically ordered, disordered and random configu- 
rations can be quantitatively analyzed with respect to 
their short-range order (SRO) correlations, and one can 
answer such questions as: which configuration is most 
similar to the random state22i22i2I or, more relevant for 
this work, which is most similar to the random site model 
(RSM)? The procedure is to calculate a lattice averaged 

spin product, n( r ^)s' f° r eacn structure (S) and cluster 
(r,i): 

]J(r, t) s = M(0)£(0) s + M(l)£(l) s + 



M(2, l)£(2,l) s + Af(8)£(8) s (3) 

Here, 0( r ^)s' is the product between multiplicities, 
M (r, t) , for r-body clusters of type t and correlation 
functions, £(r, t), for the clusters. The parameter 
r=l,2,..., is the number of sites in the cluster (l=site, 
2=pair ...8=cube), and t is an arbitrary index (e.g. 
t=2 for an r-body cluster of type 2). In a two com- 
ponent system, there is a one to one correspondence 
between symmetrically distinct clusters and (r, t) clus- 
ters (symmetrically distinct in the high-symmetry 
phase). The correlation function for any cluster is 
£(r,i) =< <Ji(jj...cr T >, where: < Oi<jj...o r > is an 
ensemble average; a = 1 for Nb; and a = — 1 for Mg. 
In general: £(0) = 1; -1 < £(1) < 1 for the site 
(point) correlation; —1 < £(2,1) < 1 for the first nn 
pair; —1 < £(2,2) < 1 for the 2nd nn pair, etc. The 
{£(r, t)} sets listed in Tabic IxTTTl arc truncated at the 
cube-approximation; i.e. only the correlations for cubes 



6 








T _ T _ T _ T _ T _ T 



FIG. 3: Filling the Mg2/ 3 Nbi /3 layers in the [111]jvt struc- 
ture with Mg (black circles) and Nb (open circles) in the di- 
rection of (a-b) [111], (c-d) [001], and (e-f) [110]. 



and their subclusters, are included. Indices in column 1 
identify symmetrically distinct subclusters of the cube in 
space group P3m, labeled as in figure El column 2 lists 
the M(r,t); columns 3-9 list the £(r,t) for the random 
state, RSM, and the [001], [111], [110], and [lll] NC c 
supercells, respectively. 

The {£(r,i)} sets for the [001], [111], [110], and 
[111] jvcc structures were calculated by counting the dis- 
tributions of distinct cube-cluster configurations in each 
structure, and the distribution of subclusters therein. 
Correlation functions for the random configuration were 
calculated using 

X ($)ijklmnop — 

= x(l)ix(l)jx(l)kx(l)ix(l) m x(l) n x(l) x(l) p , (4) 

where x(8)ijkimnop is the probability that a cube (8- 
body cluster) has configuration ijklmnop, and cc(l); is 
the probability of finding atom-z on the site labeled i in 
Figure El That is, x(l)i = 2/3 when Nb occupies site- 
i, and x(l)i = 1/3 when Mg occupies site-i. Cluster 
algebra [relationships between the x(r,t) and £(r, t)] for 
a two component system was fully described by Sanchez 
and de Fontaine—, and the extension to multicomponent 
systems was described in Sanchez et. al£&. 

The set {£(r, t)}, for the RSM, was approximated nu- 
merically by generating a 399x399x399 site simulation 
box with ideal NaCl-type ordering of Nb and Mg, ran- 



TABLE VII: The ionic coordinates (in A) for the 30-ion su- 
percells 





[OOlJjVjVM 


[IIIWt 


type 


X 


y 


z 


X 


V 


z 


Mgi 


7.98 


-0.02 


8.16 


-0.07 


7.95 


-4.06 


Mg 2 


4.00 


-0.00 


4.08 


7.95 


-0.07 


-4.07 


Nbi 


4.01 


-0.09 


8.11 


3.95 


3.96 


-4.00 


Nb 2 


3.95 


-0.04 


0.08 


3.96 


3.96 


-8.02 


Nb 3 


7.93 


-0.02 


0.00 


7.94 


3.91 


-4.08 


Nb4 


3.99 


-3.98 


4.12 


3.89 


7.95 


-4.08 


Oi 


2.04 


0.40 


8.10 


8.21 


3.90 


-6.08 


o 2 


5.96 


-0.31 


8.04 


3.82 


4.15 


-6.10 


o 3 


4.40 


1.98 


8.12 


3.91 


0.20 


-2.15 


o 4 


3.66 


-1.98 


8.10 


4.27 


7.82 


-6.08 


o 5 


2.05 


0.17 


0.00 


-0.12 


4.18 


-2.14 


o 6 


6.03 


-0.03 


12.02 


4.15 


3.84 


-2.10 


o 7 


4.18 


2.01 


12.12 


1.99 


7.85 


-4.32 


o 8 


4.00 


-1.98 


12.07 


5.95 


3.82 


-4.32 


9 


1.98 


0.00 


3.86 


5.94 


0.29 


-4.01 


O10 


6.06 


0.06 


3.96 


0.14 


10.02 


-4.05 


On 


4.02 


2.06 


3.87 


4.22 


5.96 


-4.03 


012 


4.06 


-2.02 


3.92 


3.80 


2.03 


-4.32 


O13 


0.27 


-0.01 


10.21 


2.03 


4.21 


-3.97 


O14 


4.00 


-3.92 


1.94 


-0.11 


5.94 


-4.37 


Ois 


0.11 


-0.04 


6.09 


10.01 


-0.19 


-4.32 


Ole 


3.97 


0.14 


10.11 


-2.07 


8.25 


-4.03 


O17 


4.12 


0.02 


1.95 


8.21 


1.99 


-4.00 


Ois 


3.95 


0.20 


6.13 


7.84 


-2.07 


-4.36 


Pbi 


5.75 


1.99 


10.07 


5.91 


1.70 


-6.02 


Pb 2 


5.96 


1.98 


2.41 


1.69 


5.91 


-6.03 


Pb 3 


5.84 


1.87 


6.30 


1.95 


1.92 


-1.76 


Pb 4 


1.89 


1.74 


10.10 


5.85 


5.96 


-6.11 


Pb 5 


5.95 


-2.04 


2.40 


5.85 


1.77 


-1.90 


Pb 6 


2.06 


1.83 


6.27 


1.84 


5.89 


-1.83 



domly changing 1/3 of the Mg ions to Nb, and counting 
cube-, and cube-subcluster, correlations in the simulation 

box. 

Because n( r >*)>S f° r any ordered state or random 
ensemble is a vector, one can unambiguously (except 
for truncation errors) compare average Pythagorean dis- 
tances between n( r >*)s for different states of or- 
dered. The last two rows of Table IXIIII list av- 
erage Pythagorean distances relative to the random 
configuration,i?J_ fl(mdom , and the RSM, Rs-rsm- 



Rs-rsm = ^{ [ M(0)(£(0) s 



m R SM) ] 2 



+ [M{\)(i(\) s -i(\) RSM ) f 
M(2,1)(£(2,1) s -£(2,1)rsm) ] 2 
M(2, 2)^(2,2)5-^(2,2)^) ] 2 



(5) 
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TABLE VIII: The dynamical charges of ions in [001] NC c> 



ion i 


Z* 


Z* 


7* 


ion i 




Z* 


7* 


Mgi 


2.74 


2.65 


2.65 


Mg 2 


2.74 


2.65 


2.65 


Nbi 


6.48 


6.00 


6.00 


Nb 2 


7.87 


9.11 


9.11 


Nb 3 


7.87 


9.11 


9.11 


Nb 4 


6.48 


6.00 


6.00 


Oi 


-2.61 


-3.62 


-2.82 


2 


-2.61 


-3.62 


-2.82 


3 


-2.61 


-2.82 


-3.62 


4 


-2.61 


-2.82 


-3.62 


5 


-1.99 


-7.04 


-2.09 


o 6 


-1.99 


-7.04 


-2.09 


7 


-1.99 


-2.09 


-7.04 


o 8 


-1.99 


-2.09 


-7.04 


o 9 


-2.61 


-3.62 


-2.82 


Oio 


-2.61 


-3.62 


-2.82 


On 


-2.61 


-2.82 


-3.62 


012 


-2.61 


-3.62 


-2.82 


O13 


-4.15 


-2.57 


-2.57 


O14 


-5.71 


-2.42 


-2.42 


Ol5 


-3.93 


-2.44 


-2.44 


Ole 


-5.71 


-2.42 


-2.42 


Oir 


-4.15 


-2.57 


-2.57 


Ois 


-3.93 


-2.44 


-2.44 


Pbi 


3.52 


4.41 


4.41 


Pb 2 


3.52 


4.41 


4.41 


Pb 3 


4.17 


3.93 


3.93 


Pb 4 


3.52 


4.41 


4.41 


Pb 5 


3.52 


4.41 


4.41 


Pb 6 


4.17 


3.93 


3.93 



J K 

I — I / 

Pb 



m o 



TABLE IX: The diagonal frequencies for the 30-ion supercells 
of PMN 





[001] W oc 


[111] NT 


ion 


z 


X 


V 


z 


X 


y 


Mgi 


346 


328 


340 


311 


334 


354 


Mg2 


276 


340 


344 


308 


335 


284 


IN L>1 


279 


271 


290 


289 


288 


279 


IN D2 


322 


276 


270 


272 


298 


307 


1\TU 

IN D3 


318 


286 


251 


302 


276 


269 


IN D4 


301 


288 


283 


295 


296 


291 


Ui 


385 


537 


333 


585 


271 


277 


u 2 


262 


608 


251 


675 


266 


266 




399 


331 


521 


618 


274 


256 


U4 


256 


267 


648 


548 


284 


299 


Us 


300 


640 


248 


626 


268 


247 


1 ~i 

U6 


263 


657 


257 


682 


281 


276 


U7 


297 


272 


625 


284 


670 


255 




281 


246 


633 


273 


636 


261 




292 


577 


302 


314 


587 


296 


U10 


281 


612 


278 


266 


273 


646 




285 


278 


600 


273 


268 


659 


U12 


279 


299 


613 


287 


275 


625 


Ul 3 


555 


264 


294 


281 


662 


265 


U14 


661 


267 


261 


292 


295 


639 


U15 


735 


214 


207 


290 


602 


324 


*Jl6 


557 


241 


250 


284 


628 


249 


U17 


733 


233 


253 


302 


295 


630 


Ol8 


596 


244 


245 


276 


310 


631 


Pbi 


108 


88 


76 


81 


81 


92 


Pb 2 


81 


78 


81 


87 


86 


70 


Pb 3 


93 


81 


68 


97 


76 


94 


Pb 4 


105 


76 


89 


78 


78 


81 


Pb 5 


88 


79 


70 


84 


82 


85 


Pb 6 


92 


67 


76 


86 


83 


76 



V. STRUCTURALLY STABLE STATES AND 
LATTICE DYNAMICS 



FIG. 4: 8-cite configuration. 



+ [ M(3,l)(£(3,l)s-£(3,l)H SA f) ] 2 
+ ... 

+ [A/(8)(£(8) s -£(8)«sm)] 2 } 1/2 

A smaller value of Rs-rsm implies a structure, S, that 
is more similar to the RSM, so clearly the [111]jvt is a 
better approximate of the RSM than any of the other 
supercells listed in Fig. [21 



Lattice dynamical force constants were calculated by 
the frozen phonon method; i.e. in turn, displacing each 
ion by 0.01 A in each of the three orthogonal (Carte- 
sian) directions; the force matrix is Ri a ,jf3\ where i and 
j index ionic positions, a and (3 are the Carte- 
sian directions. The dynamical matrix DM is Di a .jp = 
1/2 

Riajp/ (MiMj) , where Mi is the i-th. ionic mass. Vi- 
brational mode frequencies were found by diagonalizing 
the DM, and used to compute IR reflectivity spectra (2): 



£aP\V) = £ooa/3 + > ~n 7T 

^— ' vf - v 2 - 

i 1 
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TABLE X: Assignment of the frequencies obtained for tetragonally averaged [PMN]jvcc : n is the degeneracy, a direction 
of displacements, q the wave vector corresponding to the reduced primitive unit cell, e the averaged over the directions 
contribution to dielectric permittivity. Notice that these frequencies differ from those obtained in the straight computations 
and, hence, can be used only for a qualitative analysis. 



n 


a 


V 


q 


type 


content (e) 


3 


xyz 





(0,0,0) 


acoustic 


all ions 


4 


xy 


61 


(0, 0, 27r/3a) 


acoustic 


Pb 


1 


z 


68 


11 / n\ 

(it /a, 7r/a, 0) 


acoustic 


Pb 


2 


z 


83 


II / o /o \ 

(7r/a, ir/a, In/ia) 


acoustic 


Pb 


4 


xy 


85 


1 1 1 n In \ 

(7r/a, 7r/a, 27r/3a) 


acoustic 


Pb 


2 


xy 


94 


II I rt\ 

(ir/a,TT/a, 0) 


acoustic 


Pb 


2 


xy 


97 


(0,0,0) 


optical, Last-type 


Pb-BU6 translation (11.3) 


2 


z 


111 


(0, 0, 27r/3a) 


acoustic 


Pb 


1 


z 


121 


(0,0,0) 


optical, Last-type 


Pb-B06 translation (5.3) 


1 


z 


178 


(7r/a, 7r/a, 0) 


optical 


B, O z 


4 


xy 


192 


(0, 0, 27r/3a) 


optical 


B, 


2 


xy 


232 


/I 1 n\ 

(7r/a, 7r/a, 0) 


optical 


B, O z 


4 


xy 


249 


(0, 0, Z-K/ia) 


optical 


B, O 


2 


xy 


249 


(0,0,0) 


optical, polar 


B-O z bending (0.9) 


4 


xy 


252 


II 1 n In \ 

(n/a, n/a, Zn/ia) 


optical 


B, O 


2 


xy 


253 


(0,0,0) 


optical, polar 


a-(j X y bending (4.2) 


1 


z 


253 


(0,0,0) 


optical, nonpolar 


O xy bending 


4 


xy 


265 


(0, 0, 27r/3a) 


optical 


B, O 


1 


z 


266 


(0,0,0) 


optical 


B-Ozy bending (4.3) 


2 


z 


266 


(0, 0, Z-K/ia) 


optical 




4 


z 


270 


II 1 n /o \ 

(n/a, ir/a, Z-K/ia) 


optical 




2 


z 


284 


II 1 n /n \ 

(7r/a, 7r/a, 2n/ia) 


optical 


B, O z 


2 


z 


287 


(0, 0, 27r/3a) 


optical 


B, O 


2 


z 


292 


(k /a,K /a, 0) 


optical 




4 


xy 


292 


II 1 o /o N 

(7r/a, 7r/a, 27r/3a) 


optical 




2 


xy 


318 


(k /a,K /a, 0) 


optical 




4 


xy 


342 


(7r/a, 7r/a, 2n/ia) 


optical 


B, O z 


2 


xy 


354 


II I rt\ 

(k /a,K /a, 0) 


optical 


B, O z 


2 


z 


402 


(0, 0, 27r/3a) 


optical 


B, O 


1 


xy 


571 


(0,0,27r/3a) 


optical 


B, O 


2 


xy 


591 


(0,0,0) 


optical, polar 


B-0 stretching (0.8) 


1 


z 


599 


(k ja,K ja, 0) 


optical 


B, O z 


1 


z 


606 


(0,0,0) 


optical, polar 


B-0 stretching (0.7) 


2 


z 


660 


(k j a, k j a, 2n/ia) 


optical 


B, O z 


4 


xy 


677 


(k /a,K /a, 2n/ia) 


optical 


Ojjy 


2 


xy 


694 


(k /a,K /a, 0) 


optical 


O^y 


2 


z 


729 


(0,0,27r/3a) 


optical 


B, O 



ZLZ. 



ia i0 



= J2 Z L-y (mo/Mi f 2 



(7) 
(8) 



Here V is unit cell volume, Mi is atomic mass, m — 
la.e.m., — 5.83, the experimental electronic dielectric 
permittivity, v is the frequency of the incident radiation, 
Vi is the frequency of mode i, and a M i 7 the component 



of the DM eigen-vector for the /i-th mode, i-th ion, and 
7 direction. The damping constant was set to 60 cm -1 , 
the approximate average of damping constants that were 
fit to experimental results (Table I). 
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TABLE XI: Eigen vectors for the zone-center modes polarized 
in the z direction obtained with the help of the tetragonally 
averaged dynamical matrix of [001] jvcc 



V 


B 




O y 


o z 


Pb 


121 


0.027 


0.035 


0.035 


0.032 


-0.017 


253 





-0.072 


0.072 








266 


0.031 


-0.056 


-0.056 


-0.003 


-0.002 


606 


-0.014 


-0.019 


-0.019 


0.094 


0.000 



TABLE XII: The Wyckoff positions in the Immm group. 



a position 


ions 


4f 


Mgi, Mg 2 


2c 


Nbi 


2a 


Nb 2 


4e 


Nb 3 , Nb 4 


8m 


Oi, 3 , 4 , Os 


4i 


2 , 6 


4n 


7 , 9 , Ol4, Ol7 


4n 


Og, On, Ol2, Ol 3 


4n 


OlO, Ol5, Ol6, Ol8 


8m 


Pbi, Pb 2 , Pb 5 , Pb 6 


4j 


Pb 3 , Pb 4 



0.01 



0.00 



a>, -0.01 

>s 

CD 

L. 

CD 
c: 
LU 



-0.02 



-0.03 



0.00 0.05 0.10 0.15 0.20 
maximal Pb displacement [A] 



FIG. 5: Computed potential relief for the Last-type soft mode 
in the [001] jvcc structure. 




Mg 



FIG. 6: Atomic shifts in the equilibrium structure of [111]jvt- 
The average shift of the Pb ions is approximately along a [66c] 
direction. 
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FIG. 7: Computed pair distribution function in the 
[001]jvcc and [111]jv T structures of PMN. PDF's are com- 
pared with experimental results of Egamii— 



A. Ionic relaxation 

We first consider the relaxation of ions with respect to 
ideal perovskite positions in the Gchem structures. In ev- 
ery structure except [111]atjvm, the largest displacements 
involve Pb. In [001]^^/, the Pb-ions with 4 Mg neigh- 
bors displace by 0.47 A toward the (001)a/ s layers in 
which the O-ions are underbonded (coordinated by two 
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TABLE XIII: Multiplicities and correlation functions for: a random distribution; the random-site model (RSM); the five ordered 
supercells. 



Cluster 


Mir tV 


Random 


RSM t 


rooii 


fllll 
[iiij 


[1101 

[11UJ 


fllll fjnr'i 
[i i ij jVCC 


fllll WT 
[i i ij iv 1 


0* 


I 


I 


i onofin 


I 


I 


I 




I 


1 


1 
1 


-1/0 


-u.oooou 


-1/0 


-1/0 


1 /*} 

-1/0 


1 /■} 


1 /*} 
-I/O 


!J 


3 


1 /9 
i/y 


\J . 1UUOU 


5/9 


-1 /3 
1/ o 


1 /Q 
1/ J 


-1 /3 

1/ o 


-1 /q 
i/ j 


ik 


u 


1 /Q 
i/y 


u.uzyou 


1 /9 
i/y 


1 /3 
1/ 


i /q 
-1/ J 


/Q 

U/ C7 


n 


lO 


A 


1 /Q 
i/y 


u.uuy± £ i: 


1 /■} 

-1/0 


n 


1 /^ 
1/0 


1 /■} 
-1/0 


1 /^ 

1/0 


iilf 
!J K 


1 9 


1 /97 


u.uoiyu 


1 /■} 
-1/0 


1 /■} 
1/0 


1 /Q 

-i/y 


1 /q 
i/y 


9 /q 
z/ y 


ikn 


Q 
O 


1 /97 

-1 / Z 1 


n n^i no 


1 /■} 


1 /■} 
-1/0 


i /*} 
i/o 


1 /■} 
-1/0 


n 


ijo 


z4 


-1/2Y 


-U.UU044 


-1/3 


U 


i /n 
1/9 


i /a 
l/y 


U 


1JK1 


Q 
O 


1 /SI 

1/ Ol 




1 

1 


1 /■} 
-1/6 


/q 
o J y 


1 /q 
i/y 


1 /q 

-i/y 


ikmo 


1 9 

1Z 


1 /SI 


n n p i^9s 

U.Ui_MZiO 


1 

1 


1 /■} 
1/0 


1 /q 
l/y 


1 /q 
i/y 


1 /q 
-i/y 


ijlm 


Q 

o 


1 /Q1 

l/ol 


-U.U4oyo 


1 /■? 
-1/3 


1 /■? 
-1/3 


i /Q 
1/3 


1 /■? 
-1/3 


U 


iknp 


;[ 


1 /SI 

1/ Ol 
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1.080 


1.201 
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* M(r,t) is the multiplicity, per site, of the r-body cluster (r=l,2,...8 sites) of type t. 

* Real numbers are from a 399x399x399 site simulation, and fractions are inferred, from multiple simulations. 

* The Zero- or empty cluster. 

* Rs-Random is the average (divided by the number of correlation functions, 22) Pythagorean distance between the LASP of 
structure S, and the LASP for a random alloy; and similarly for Rs-rsm- 



Mg 2 + ionai^iS). In [1W] NNM , the Pb-ions with 2 Mg 
neighbors displace by 0.47 A toward the underbonded 
O-ion between the 2 Mg. The other structures do not 
contain nearest neighbor Mg-Mg pairs; nonetheless, each 
Pb-ions in an asymmetric environment of Mg and Nb 
is displaced off-center by a nonzero "local field" i*^> The 
magnitude of these displacements ranges from 0.25 A to 
0.27 A for the [001]atcc" an d [Hl].iVT structures, but is 
only 0.04 A for the [111] nnm structure. The small Pb 
displacements in the [111]awm structure result because 
of a special local environment of Pb in this structure. 
The largest ionic motion is by the Nb-ions, which shift 
0.09 A towards nearest neighbor (lll)M g planes. In the 
[111]att structure, Pb deviates from the [111] axis; a net 
displacement of (0.17,0.17,0.08) A occurs. 

Next, we consider the nature and result of symmetry- 
breaking relaxations leading to the Gsss structures. All 
the 15-ion supercell have multiple instabilities. While 



the different shapes of the supercells studied prevents 
direct comparison of the modes between the different 
cells, there are both ferroelectric instabilities of the Last 
type (Pb motion against the other ions), and antifer- 
rodistortive (tilting) instabilities dominated by O mo- 
tion. Fig. shows the potential barrier for the Last 
mode in [001]atcc" structure. With respect to the primi- 
tive perovskite cell, every supercell shows wide dispersion 
of Pb- and O- dominated instabilities across the Brillouin 
zone. Nb motion opposite O in Nb-O-Nb... chains and 
planes also plays a role in the FE instabilities. We note, 
however, that long Nb-O-Nb... chains and planes are un- 
likely to form experimentally. 

Given the existence of multiple instabilities in PMN 
supercells, it is not surprising that the Gsss states in- 
volve the freezing of multiple modes. Figure shows the 
complexity of the local structure in fully relaxed [111]jvt- 
Pair distribution functions (PDF) for the [001]jvcc m 
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PA/nmm (dotted) and PI (solid) structures are shown 
in Fig. [7| The peak at 2.8 A is already split by the differ- 
ent Pb-displacements in the Gchem structures, and the 
splitting is significantly enhanced in the Gsss structures. 
The agreement between the predicted and experimental 
PDF is quite good. 

We estimate the polarizations of our stable structures 
by using the displacements of each ion from its posi- 
tion in Gchemi and estimated Born effective charges de- 
termined from averaging the calculated Born effective 
charges for the [001] jvcc structure 



38,39 . 



J Pb 



2.6; Zi 



7.4; Z? 



4.0; 
-2.5. 



^Mg ~ ^Nb ~ >-^> ^Oy ' Ox 

The results are shown in Table IVAl The estimated po- 
larization magnitudes are similar in magnitude, though 
about 10 % smaller than the zero-temperature polariza- 
tions estimated from first principles for PbTiOs— and 
PbSc!/ 2 Nbi/2 03 4i Remarkably, three of the five struc- 
tures have SSS polarizations that are pseudomonoclinic 
(close to an [aab]-type direction). 



TABLE XIV: Estimated polarization (in C/m 2 ) of struc- 
turally stable states of PMN supercells. 



structure 


P. 


Py 


P z 


|P| 


[001]jvjvm 


0.525 


0.191 


0.000 


0.558 


[1 10] JV JVM 


-0.158 


-0.158 


-0.547 


0.591 


[111] JV JVM 


0.474 


0.476 


0.016 


0.672 


[001] NGC ' 


-0.228 


-0.243 


0.385 


0.509 


[111] jvt 


-0.238 


-0.221 


0.346 


0.474 



B. IR spectra and lattice dynamics of relaxed 
structures 

The IR reflectivity spectra are shown in Fig. |H1 and 
compared with our experiment. It is seen that the best 
agreement between the theory and experiment is for 
[111] jvt structure. The complex dielectric spectra in 
the phonon frequency range are plotted in Fig. to- 
gether with experimental data obtained. We have qual- 
itative agreement for high frequencies. At low frequen- 
cies, experiment shows very high values of permittivity 
connected with relaxations of PNC. These frequencies 
are below phonon's. Relaxation contributions to dielec- 
tric permittivity were not included into our theoretical 
consideration (we discussed these problems in a previ- 
ous publication 42 ' 43 ). Notice different scales in the ex- 
perimental and theoretical plots. The total zone-center 
phonon density of states of the [lll]jvcc" an d [111]jvt 
structures are shown in Fig. 1101 It is seen from this plot 
that the high energy phonon band is not very dispersive 
compared to the low frequency bands and compared to 
the bands known for ABO3 perovskites. This can be ex- 
plained by folding the bands due to ordering of Mg and 
Nb. 



TABLE XV: The ionic coordinates for the structurally stable 
[111] jvt supercells. In units of the basic vectors. 
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Mode assignments were then made at several levels 
of approximation by averaging the [111] jvt DM. 4 — The 
lowest-frequency modes of the supercell are related to the 
folded acoustic modes of the parent 5-atom cell: exclud- 
ing zero frequency zone center acoustic modes, eigenvec- 
tors of the low-frequency modes in the frequency interval 
19-100 cm -1 are significant only on Pb. This corre- 
sponds to extremely low-frequency acoustic antiphase Pb 
displacements, which arise due to a very small diagonal 
element of DM for Pb. 

Next, to understand how local 1:1 order may affect 
the phonon properties, the [111] jvt dynamical ma- 
trix was averaged with respect to a 1:1 ordered struc- 
ture. The [111] nt structure is the only one studied 
which is commensurate with such averaging. Effectively, 



the [111]m S2 



3 Nb x 



/3 



layers were treated as though they 



contained only one average T-cation, T=(Mg 2 /3Nb 1 / 3 ). 
Symmetry analysis of [111]jvt in Fm3m, yields the fol- 
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Wavenumber (cm 1 ) 



FIG. 8: Computed (a-e) and experimental (f-g) IR reflectiv- 
ity of PMN: a) [Q01] NNM ; b) [0U] NNM ; c) [111] NNM ; d) 
[Q01}ncc>; (e) [H1]jvt; f) 20 K; g) 300 K. Dotted lines are 
damped oscillator model fits (see Table HJ to the experimen- 
tal IR reflectivity spectra (solid lines) and THz spectra (solid 
points) . 



lowing modes: A lg (R) + E g {R) + F lg (s) + 2F 2g {R) + 
4Fi u (IR) + F2u(s). Here, IR and R indicate IR- and 
Raman- active modes respectively, and s indicates silent. 
Computed eigenvectors of the [111]jvt -PI structure were 
projected onto vectors of the averaged computation to 
obtain the modes expected to have the highest infrared 
and Raman activity. The results are shown in Fig. II 11 
and compared with experimental data from inelastic neu- 
tron scattering^ Raman scattering and FIR reflectivity 
measurements (this study). 

Polar TO modes [F lu (IR)] are labeled TOl, T02, 
T03, and T04. The relative contributions from Pb, T, 
Nb, 0| h and 0_l in the TOl mode (100 cm" 1 ) are: -0.16, 
0.21,0.27, 0.32, and 0.41 A; where: O M is displaced along 
B-0 bonds; Oj_ is displaced perpendicular to B-0 bonds. 

The T02 mode at 270 cm" 1 , is predominately of the 
Slater type. Neutron^ and our IR measurements (Ta- 
ble I) exhibit bands near this frequency. The T03 mode 
near 340 cm" 1 is from Nb moving opposite to T. Ex- 
perimental evidences for this band is not definitive, al- 
though a mode is found near this frequency in the fit 
to room-temperature IR data. The T04 mode near 580 
cm" 1 originates from 0||-(T,Nb) stretching, and it is 




Frequency [cm" 1 ] 



FIG. 9: Experimental (a-b) and computed (c-d) dielectric 
permittivity and losses: (a-b) squares denote experimental 
microwave measurements obtained by using a THz spectrom- 
eter, circles IR measurements; the dashed and solid lines are 
the results of the fits to IR reflectivity at 20 K and 300 K 
respectively; (c-d) the dashed and solid lines correspond to 
the [001]jvcc and [111]jvt structures respectively. 

clearly observed in the IR and neutron data. 

In first-order Raman scattering from Fm3m crystals, 
A lg and E g modes should be visible in VV geometry, 
while F-ig modes should be visible in VH geometry. The 
calculated A\ g Raman active mode at 810 cm" 1 matches 
quite well the intense peak that is observed experimen- 
tally at 780 cm" 1 in VV spectra. In [111]att Fm3m, this 
mode corresponds to a fully symmetrical, IR inactive, O- 
breathing mode. The calculated Raman active E g asym- 
metrical breathing mode near 590 cm -1 is broadened due 
to local symmetry breaking, and corresponds with exper- 
imentally observed bands at 500 and 600 cm -1 . In all 
cases, we had some frequency intervals (and even a mul- 
tipeak structures) when projecting the modes obtained 
for GSS onto the modes derived from symmetrized DM. 

The calculated Raman active -F2g(2) mode near 360 
cm" 1 involves mostly symmetrical Oj_ displacements. A 
wide envelope of the bands is seen in our Raman exper- 
iment at this frequency. The calculated low-frequency 
Fi g Raman mode near 80 cm" 1 is mostly from Pb- 
displacements. In the VH spectrum, the low-frequency 
band is a doublet (45 and 62 cm" 1 ). Both compo- 
nents are close to the zone-boundary frequency of the 
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200 400 600 800 

frequency [cm 1 ] 



FIG. 10: Computed zone-center phonon density of states in 
structurally stable states. 



TA branch (50 cm 1 ) according to neutron scattering 
data££ (Fig. [TTJ>- 

Experimentally, there are strong Raman peaks near 
150 cm -1 and 250 cm -1 . These peaks probably arise 
from local symmetry breaking (octahedral tilting and/or 
polar distortions), and not from Nb-T ordering. 



VI. DISCUSSION 

All the PMN supercells studied here (Fig. |5J) are 
unstable with respect to lower energy structurally sta- 
ble states, i.e AEchem > AE$ss- Relaxation ener- 
gies, AEsss — AEchem, are large relative to the 
thermal vibrational energy, which suggests that local 
displacive symmetry breaking should persist to high 
temperatures, generating local strain fields. Some in- 
stabilities, however, have barriers of the order of ther- 
mal energy, especially Last-type modes corresponding to 
Pb-BC>6 stretching vibrations in finite-size chemically 
ordered regions. The softest vibrations are exhibited by 
Pb-ions that are coordinated by a highly symmetric array 
of Mg- and Nb-ions, or Pb-ions surrounded by Nb-ions 
only. Such Pb-ions have diagonal Pb-frequencies (the ra- 
tio of the spring constant and mass) of only ~ 80 cm" 1 . 

Of the supercells studied here, [111]jvti m 
Gsss = -PI, has the lowest energy, and with re- 



spect to chemical SRO, it is the best approximate of the 
RSM. This suggests that it is also a good approximate 
for chemically ordered domains in a PMN crystal with 
PrriSm global symmetry. EgamiaS and Naberezhnov 
et al£Z> reported that the local ferroelectric instability 
in PMN below Burns temperature is associated with 
Pb-displacements in opposition to the other ions and 
this is essentially what occurs in the [111]att Last-mode 
dominated FE instability, Immm — > PI. Note, 
however, that Last-type instabilities occur in all five 
supercells. Computed distributions of the interatomic 
distances in the 30-ion supercells show that Pb-0 
distances split into two main groups at 2.5 A and 3.2 
A , in agreement with experimental 

Computed dynamical charges in PMN are not as large 
as is typical for simple ABO3 perovskites 39 : the Nb 
charge varies from 6.0 to 9.1. An FE-active ion that is 
surrounded by less FE-active ions, e.g. Nb surrounded by 
Mg, typically has a reduced dynamical charge, e.g. rela- 
tive to Nb surrounded by Nb; and the dynamical charge 
of the less FE-active ion increases (cfi££) . This happens 
because charge transfer is reduced when an ion of rela- 
tively low electronegativity is surrounded by ions with 
higher electronegativities 4£ 

Computed phonons for all the structures shown in 
Fig. IHlexhibit three similar features: (1) B-O-B stretch- 
ing modes at 500-900 cm" 1 ; (2) mixed B-O-B bending 
and O-B-0 stretching modes at 150-500 cm" 1 ; (3) Pb- 
B06 stretching modes at v < 150 cm" 1 . Differences 
between the computed IR spectra for different super- 
cells (Fig. |8J) are purely quantitative, not qualitative. 
The [111]jvt spectrum is most similar to the experimen- 
tal one: the 500-900 cm" 1 band is relatively narrow, 
as in the experiment (Fig. |HJ); but the other computed 
bands are less intense than their experimental counter- 
parts. This is an artifact of using the same damping 
constant for all computed modes (60 cm -1 is an aver- 
age of the experimental values). The 60 cm" 1 damping 
constant, which we used in our computation, is signifi- 
cantly larger than the usual value for pure perovskites, 
^20 cm" 1 . A large value for disordered PMN is expected 
because the random fields caused by chemical disorder 
broaden the phonon density of states. Experimentally, 
however (c.f. Sec. some modes exhibited extremely 
small damping constants, which implies sharper features 
in the reflectivity spectrum. 

One can see in Fig. 1111 that the number and distribu- 
tion of the modes in Immm group over frequencies are 
much richer than in FmSm geometry. For example, the 
single A\ g peak in Fm3m group can be connected with 
three A g modes in Immm group (857, 800, and 759 
cm" 1 ), in agreement with Raman spectra. The distribu- 
tion of the frequencies of these peaks in Immm group 
is due to differences in frequencies for oxygen ions in dif- 
ferent Wyckoff positions: the stretching vibrations of the 
oxygen ions in Nb-O-Nb- chains are harder than other 
oxygen stretching vibrations, due to ferroelectric insta- 
bility in these chains (combined with nearest Pb ions) 
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FIG. 11: Comparison of low-temperature experimental data and computational results: a) 12 K neutron scattering data™ (full 
triangles and squares correspond to the longitudinal and transverse phonon dispersion branches; the bands show the frequencies 
of continual neutron scattering); b) Raman scattering data in two geometries (VV, solid line, and VH, dotted line) recorder at 
77 K; c) dielectric loss e"(y) obtained from the fit of IR and THz spectra at 20 K; d) Results of first principles computations 
in Fm3m symmetry; e) Mode assignments for PMN [111]att structure in space group Immm: two last columns correspond 
to different polarization of Raman spectra. 



and a large shift of oxygen ions in the equilibrium state 
(see Table Hxjl . There is also a difference between ionic 
vibrations in Mg-O-Nb bridges and in Nb404 squares in 
xy plane. 

Ideally, A g modes should be seen only in VV geom- 
etry. The presence of some intensity in VH geometry 
can be regarded to local fields and distortions. Notice 
the existence of an IR active mode (704 cm" 1 ) in the 
same frequency interval which is due to a large (mini- 
mal) size of the supercell in comparison with the unit 
cell at FmZm symmetry. 

At the frequency of E g mode of FmZm group, 
there are several Raman active symmetrical breathing 
A g (663, 601, and 534 cm -1 ) and one antisymmetri- 
cal breathing (599 cm -1 ) modes. The splitting and the 
absence of large intensity in VH geometry in experimen- 



tally measured spectra in this frequency interval testifies 
against E g mode in Fm3m group but agrees quite well 
with the computed splitting in Immm group and with 
A g symmetry of three of the four modes. In the same 
frequency interval there is a wide distribution of IR active 
modes that well corresponds to the wide peak in the FIR 
spectrum (T04 mode in terms of Fm&m symmetry). 

The Fm3m group assignment can hardly explain the 
presence of a peak in Raman and IR spectra at a fre- 
quency a little bit higher than 400 cm -1 . The compu- 
tation within Immm symmetry does give both Raman 
active (A g , 443, and B 2g , 423 cm -1 ) modes as well as IR 
active modes (412, 418 cm -1 ). 

The position of F 2g (2) mode in Fm3m group is no- 
ticeably higher than the position of the intensive peak in 
Raman spectrum in VH geometry. The Immm group 
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gives a few suitable modes in the right place and gives 
a wide distribution of these modes over frequencies that 
corresponds to a large width of the Raman lines observed. 
However, at the present time, we cannot prove that the 
peak of the Raman intensity in the computation corre- 
sponds to the peak in the experiment, because we did 
not compute the intensities. 

Finally, we want to notice that there are a few low 
frequency modes spread over the frequency interval 25- 
96 cm -1 (below the frequency of the soft mode) which are 
Raman and IR active. These modes appear because of 
comparatively low symmetry of the supercell. All they 
are in the range of the acoustic modes in the reduced 
(5-ions) unit cell and are connected with antiphase Pb 
displacements. 

Taken together, the first principles results, and ex- 
perimental data from neutron scattering^ 45 Raman, and 
IR reflectivity enables assignments of specific vibrational 
modes to experimental phonon peaks. Following Akbas 
and Davies^ PMN is approximated as the RSM. Posi- 
tions of IR-active, Raman-active and silent modes, were 
calculated, and distribution functions for all these modes 
have finite width, and some are split. 

There are a few low frequency modes in the range 25- 
96 cm -1 (below the frequency of the soft mode) which 
are both Raman and IR active. They reflect antiphase 
Pb-displacements 4 ^ that are allowed by the low supercell 
symmetry, and are all in the acoustic-mode range for the 
reduced (5-ion) unit cell. 

VII. CONCLUSIONS 

We studied five ordered supercells of PMN. Instabili- 
ties involving Pb off-centering are observed in all super- 



cells and instabilities involving the tilting of O octahe- 
dra centered on Mg ions are observed in most. When 
fully relaxed, the [111]jvt structure described here has 
the lowest known first-principles energy for an ordered 
PMN structure. It is commensurate with the random 
site model (RSM) of PMN 8 , and has short -range Mg-Nb 
correlations that are most comparable to the RSM. Its 
Raman and IR spectra and pair distribution functions are 
in qualitative agreement with experiment, allowing mode 
assignment of the dominant features of the spectra. In 
its lowest-energy state, it is polarized in approximately 
along an [bbc] type direction, yielding pseudomonoclinic 
symmetry. Unphysical aspects of this perfectly ordered 
[111] att supercell as an approximant for the 1:1 ordered 
regions in experimental PMN, include the lack of a dis- 
ordered matrix, the fact that there is no experimental 
evidence for ordering on the Mg 2 /3Nb!/ 3 sublattice, and 
the presence of infinite Nb-O-Nb... chains (which influ- 
ence the lattice dynamics). 
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